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Abstract 

With the help of the Mellin-Barnes transform, we show how to simultaneously re- 
sum the expansion of a heavy-quark correlator around = (low-energy), = 4m^ 
(threshold, where m is the quark mass) and q^ — > — oo (high-energy) in a systematic way. 
We exemplify the method for the perturbative vector correlator at ©(a^) and ©(a^). 
We show that the coefficients, ri(r7,), of the Taylor expansion of the vacuum polarization 
function in terms of the conformal variable uj admit, for large n, an expansion in powers 
of 1/n (up to logarithms of n) that we can calculate exactly. This large-n expansion 
has a sign-alternating component given by the logarithms of the OPE, and a fixed-sign 
component given by the logarithms of the threshold expansion in the external momentum 



1 Introduction 



Two-point correlators of heavy quarks are very useful objects for extracting parameters 
of QCD such as, e.g., the quark masses and the strong coupling constant a^, both from 
experiment [11 El E] and also from the lattice[4]. However, even our perturbative knowledge of 
these functions is quite limited due to the complexity of the required diagrammatic calcula- 
tion. For instance, in the case of the vector vacuum polarization function, only at 0{ag) and 
0{al) is the analytic result fully known [5j. At ©(a^), state-of-the-art techniques produce 
partial answers for this function only in the form of local expansions around three particu- 
lar values of the external momentum, to wit = 0, (7^ — >• —00 and = 4m^ where m is 
the heavy-quark mass. These correspond to the so-called low-energy expansion, the operator 
product expansion and the threshold expansion of the correlator. 

A phenomenal effort has been devoted to the calculation of an ever-increasing number of 
coefficients in all of these expansions. For instance, for the 0(0^) contribution to the vacuum 
polarization, no less than 30 terms are known from the low-energy expansion [6(i7j; although 
in most of all other situations this number is reduced to just a few [8ll9lll0 1 [TT | ll2 | [T3]. Ideally, 
of course, one would like to be able to know the full function (including its imaginary part) 
and, consequently, a strategy has to be devised in order to reconstruct the function from the 
three expansions, even if this is only approximate. 

Up to now, the method that has been employed in this reconstruction is that of Fade 
Approximants |14) . Fade Approximants are ratios of two polynomials whose coefficients 
are tuned so that they reproduce the Taylor expansion in around the origin of a given 
correlator, to the highest possible degree. However, even though the vacuum polarization 
is an analytic function and, therefore, admits a Taylor series around q"^ = 0, this is not so 
either at threshold or for — )• — 00 due to the presence of non- analytic dependence in the 
form of square-roots and logarithms of the momentum. Consequently, all by themselves. 
Fade Approximants cannot approximate the correlator in the full complex plane, unless this 
nonanalytic dependence is somehow removed. In order to do this, a battery of judiciously 
chosen functions is used for subtracting the nonanalytic dependence on q"^ from the original 
correlator. Fade Approximants are then constructed for the subtracted correlator, rather 
than for the correlator itself. As pointed out in, e.g., ref. |131 115) this procedure has a certain 
degree of arbitrariness. This ambiguity in the procedure can potentially be the source of a 
systematic uncertainty, and it is important that this uncertainty be quantified. Normally, 
this is attempted by varying among several of the choices made in the aforementioned Fade 
construction. The results obtained with the Fade method are certainly interesting but, as 
higher and higher precision is always in demand, we think it is also interesting that one could 
devise alternative methods which could confirm (or if need be, modify) these results. 

Therefore, we think it would be desirable to have an alternative procedure which could 
reconstruct the full function from the three expansions in an unambiguous way, and with 
which we could have a better handle on the possible systematic uncertainties involved. In 
this paper we will present what we believe to be one such method. Furthermore, unlike Fades 
which are essentially numerical, our reconstruction is totally analytic. 

The main observation is that, in the conformal variable u in which the q'^ cut complex plane 
is mapped onto a unit disk, the coefficients of the threshold expansion and the OFE determine 
the asymptotic behavior of the coefficients of the low-energy expansion. This property is 
easiest to recognize in the language of Mellin transforms, thanks to a mathematical result 
which is known as the Converse Mapping Theorem [16j. 

An intuitive way to understand this is the following. A finite set of Taylor coefficients 
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knows nothing about the region of analyticity of the function. So, for example, even if you 
know a thousand Taylor coefficients, strictly speaking, you still do not know the value for the 
radius of convergence. The reason is because the radius of convergence is a limit procedure and 
things can change completely from the coefficient a thousand and 1 onwards. Since the region 
of analyticity knows about the asymptotic tail of the Taylor coefficients, it is not strange that 
by giving precise information on the properties of the function at the boundary of this region 
of analyticity (i.e. the threshold expansion and the OPE) one may gather information on 
this asymptotic tail of the Taylor coefficients. How to make this statement useful beyond just 
words, is what we have tried to do in this work. 

As it turns out, the asymptotic behavior of the low-energy coefficients sets in so quickly 
that already for the second or the third coefficient it is a good approximation. This results 
in an efficient way to compute them. Once all the Taylor coefficients are known (even if only 
approximately) one can resum them and reconstruct the full function. As we will see, this 
resummation expresses the result for the full function as a unique combination of polylog- 
arithms (and derivatives thereof) plus a known polynomial. Although we have set up the 
problem within the context of a two-point correlator of a heavy quark, we think that the 
method we have developed can be of more general use, and could be applied also in other 
contexts where several expansions to a function with a given analytic structure are known. 
Mellin transforms have also proven very useful in other contexts in perturbative calculations 



In the next section we will set up the necessary definitions and describe in detail our 
method and the relevant mathematical results it is based on. In section 3, we will test the 
method with the case of the 0{a1) vacuum polarization function which, although not fully 
known analytically, is sufficiently known to allow interesting nontrivial checks. In particular, 
since 30 low-energy coefficients have been calculated up to now [H [7], we will be able to 
compare our prediction for them with their exact value. In section 4, we will apply the method 
to the 0{ag) vacuum polarization function, for which only a few terms in the low-energy, 
threshold expansion and OPE have been calculated (H U\ [HI [El [21]. We will reconstruct 
the full function, make again predictions for the coefficients of its low-energy expansion, and 
compare them with those already present in the literature. The final section will be devoted 
to some conclusions and outlook. Finally, in the appendix, we compare our approximation to 
the exact result in the case of a known function, choosing the result at 0{a^) as an illustrative 
example. 

2 Description of the method 

Let Il{z) be a function of the complex variable z with the following properties: 

• It is analytic in the disc |2;| < 1 where, consequently, it admits a convergent expansion 
when z — )• as 



which, from now on, we will call the "Taylor" expansion. 
• It has a cut along 1 < Re z < oo. In the limit z — )• 1, it can be expanded in the form 



[I3ll8l[l9l|2n|. 



oo 




(1) 
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n(z) ~ Yl ^(P^ ^) (1 - - ^) ' (2) 

where A; are integers but p may be integers or half-integers. This expansion wiU be called 
the "threshold" expansion. 

• Finally, it admits an expansion for z — )• — oo such as 

n(z) ~ Vi?(p,A;) llog'^Mz) , (3) 

2— — ^ ' 

p,k 

where p and k are both integers. We will refer to this expansion as the "OPE". 

The problem we are trying to solve is how to reconstruct the function n(z) which matches 
onto the above three expansions in the corresponding regions, in an analytic and systematic 
way. That this is possible is already guaranteed by the Taylor expansion (which is conver- 
gent) and the theorem on analytic continuation. However, an exact reconstruction requires 
complete knowledge of the infinite set of coefficients C{n). Realistically speaking, this type 
of information is certainly not available. When only a finite number of C(n)'s is known, one 
clearly has no other option but to resort to the OPE and Threshold expansions to try to 
reconstruct the original function 11(2:). It seems reasonable to imagine that the more terms 
one knows, the better one will be able to determine this function. However, to the best of our 
knowledge, a concrete method to analytically reconstruct n(z) is not available in the litera- 
ture. Therefore, we would like to explain how one can do this reconstruction in an analytic 
and systematic way. This we will do with the help of the Mellin transform, using the vector 
two-point correlator for a heavy quark as an example, and illustrating how the method works 
in practice. 

Let us recall the definition of the vacuum polarization function n(g^) through the corre- 
lator of two electromagnetic currents j^{x) = q{x)'y^q{x), where q{x) is a heavy quark, 

{g,uq' - q,q.) Ho") = -ijd^x e'^^ { |T j^(x)i,(0)| ) , (4) 

and where q'^ is the external four-momentum. In QCD perturbation theory, n(g^) may be 
decomposed to 0{ag) as 

n(g2) =nW(g2) + (^) nW(g^) + (^)'n(2)(g^) + (^)' n(3)(,2) +o(af). (5) 

For definiteness, Og denotes the strong coupling constant in the MS scheme at the scale fi = 
n^-pole, but the precise definition is not important for the discussion which follows. Equation 
([5|) will be understood in the on-shell normalization scheme where a subtraction at zero 
momentum has been made in such a way as to guarantee that 11(0) = 0. As it is well known, 
the vacuum polarization in Eq. ([5]) satisfies a once subtracted dispersion relation, i.e. 

ntf) = ^7„ -imn(t + fe). (6) 

When massless cuts are disregarded (as we will assume henceforth), the spectral function in 
([6]) starts at t = Am?, where m is the heavy quark pole mass. In this case, the function n(g^) 
may be expanded for — t- in powers of 
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Figure 1: Conformal mapping between z and m, Eq. ([8]). 

and it admits a Taylor expansion, exactly as in Eq. ([1]). The other two expansions, Eqs. (j2|3p 
coincide with the threshold expansion (for — )• 4m^) and the OPE (when — )• — oo). 

In order to bring the three expansions ([UlS]) to a common ground, it is very convenient to 
make the conformal transformation 



4a; 



(l+o;) 



l-VT 



i + \/r 



(8) 



which maps the cut z plane into a unit disc in the to plane, as we can see on figure [TJ The 
cut z £ [l,oo[ is transformed into the circle |a;| = 1 and the points z = into w = 0, z = 1 
into CO = 1, and the limit z — )• +oo it ie into uj — )• —1 it ie, and z — )• — oo into uj — )• —1. 

In terms of the new variable u, the vacuum polarization function n(a;) also obeys the first 
and second expansions ([HEl), while the OPE becomes formally like the threshold expansion, 
except for a flip of sign in u. To wit: 



The Taylor expansion becomes 



n) UJ 



(9) 



n=0 



whereas the OPE (w — )• —1) and threshold (lo — )• +1) expansion become 



n(w) 



^ J]l^W(A,p) (iT^)MogP(lT 

X,p 



(10) 



The two sets of coefHcients, fl^~^\X,p) and fl^~\\,p) are, in principle, unrelated as they 
correspond to expansions of the same function at two different points. However, one can find 
the connection between the two sets of coefficients C(n) in Eq. ([TJ and 0(n) in Eq. Q. They 
read 
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CM - ^("^^(^ + ") V ^(^)^ (12) 

As we will now see, the coefficients Q^^\X,p) in (|10p determine the behavior of Q{n) for 
large values of n. In order to see this, it is very convenient to redefine the conformal variable 
u as 

uj = ±e"* . (13) 

The limit t — )• with the "+" sign corresponds to the threshold expansion, whereas the "— " 
sign corresponds to the OPE. Then, one can make use of the integral representation 

c+joo 

T{s) = / dt f-^ e"* e"* = / — t-' r{s) , (14) 

Jo J 2?7r 

c— ioo 

where the integral on the righthand side term runs over a straight line in the s complex plane, 
with c any constant real number, < c < cxd. The second identity allows one to express the 
vacuum polarization function ([9]) as an inverse Mellin transform in the t variable as 

c+ioo ^ 
f' li 

n(±e-*) =^ / ^ r{s) J2 (±1)" ^i^) ■ (15) 

c-ioo "=1 

In this expression the integration path is a vertical line in the s complex plane located at 
Re{s) = c, within the region of full analyticity of the integrand. This region is usually called 
"the fundamental strip". 

The usefulness of the expression (jlSp for our purposes is due to a relatively recent math- 
ematical theorem, known as the Converse Mapping Theorem |16) . which assures that the 
t — )• behavior of the function n(ibe~*) is completely determined by the singularities in the 
s complex plane of the function appearing in the integrand of (llSp . i.e. 



r(s) G±(s) where G±(s) = ^ (±1)" 0(n) n"^ . (16) 

n=X 

Furthermore, these singularities are just limited to isolated poles, i.e. the above combination 
r(s) G±{s) is a meromorphic function. In general, the Converse Mapping Theorem establishes 
the following dictionary [li3j between a given function f{t) and its Mellin transform Ai[f]{s): 

c+ioo 

M[f]{s) = dtf-' fit) ^ f{t)= j —t-^M[f]{s) (17) 

c—ioo 

p,fc p,k ^ ' 
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In the expression above the symbol x means "singular expansion" [16], i.e. the sum of all the 
negative powers of the Laurent expansions of the function around every pole 0. In particular, 
we would like to call the reader's attention to the result shown in Eq. (jlSp . which says that 
the expansion of the function at t — )• is given by the pole p (together with its multiplicity 
k), and the residue Tp^^, of its Mellin transform. This property will allow us, for instance, to 
compute the behavior of the coefficients Q(n) in the expansion ([9]) for n large. 

One must first realize that the position and residue of a pole of the G±{s) function is 
entirely given by the asymptotic behavior as n — )• oo of the Q{n) coefficients. This is very 
important since it means that the function Q{n) admits an expansion in powers of 1/n for n 
large. 

For instance, if 

0(n) = -3^1og'=n, (19) 

one finds that 

G+is) = i-l)'^ C^'\s + p) , G.is) = {-l)^+'r^^'\s+p) (20) 

where (^^'^'^ (s) is the A:-th derivative of the Riemann C function, and r/^'^) (s) is the k-th derivative 
of the Dirichlet r] function, 

oo oo 

C{s)=J2^" and r,{s) = -J2i-^r n-'= {1-2^-^) C{s). (21) 
n=l n=l 

Now, 



C(^)(.)x(-i; 



. nk+i) 



{s - 1)^^+1 ' 

exhibiting the pole at s = 1. So the only pole of ^4.(5) in Eq. (I20p is located at s = 1 — p. 
One immediately sees that the location of this pole is correlated with the power of n in the 
large-n behavior of the J7(n) coefficient (llOp . while the multiplicity of the pole is given by the 
power of the log n. On the contrary, the Dirichlet r] function has no pole. Of course the results 
in G+{s) and G_(s) get interchanged if the i}{n) coefficient has the same large-n behavior 
but with an alternating sign. 

Looking at the vacuum polarization function 11 in Eq. (llSp . one sees that one goes from 
the threshold expansion to the OPE by the addition of the alternating sign in the definition 
of the function G±{s). This means that both expansions will be simultaneously encoded in a 
0,{n) coefficient admitting an expansion for large n of the form 

1 k 

0(n) ~ J7^^^(n) = ^[ap,, + (-l)«/3p,,]^. (22) 

p,k 

In fact, the coefficients ap^k encode the information from the non-analytic terms of the thresh- 
old expansion, and the /3p^fc the information from the non-analytic terms in the OPE, as we 
will see explicitly in the next sections. 

To illustrate the procedure, we will ffist treat the case of the Il^'^\q'^) contribution to 
the vacuum polarization, Eq. ([5]). Since no less than 30 terms are known in the Taylor 
expansion of this function ([1]), this information will be very helpful to check our results and 
our understanding of the systematic errors involved. Furthermore, we will approximately 



^Sometimes this expansion is also called the "principal part" of a meromorphic function. 
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reconstruct the function in the complex plane, including its imaginary part, and will be able 
to compute the value of a constant in the threshold expansion of this function which, up to 
now, has not been possible to obtain by ordinary diagrammatic methods. After that, we will 
move to the more interesting (and also more difficult) case of Il^^\q'^), where we will repeat 
the same analysis, except that now only 3 coefficients from its Taylor expansion ([1]) together 
with just a few from the OPE and threshold expansion, are currently known. 

3 The case of n(2)(g2) 

Let us begin by considering the function H^'^^ . The threshold expansion ([2]) has been calculated 
in [121 [15] to be 



n(2)(z)^^^^i^=M + |A(0,2) log2(l-z) + A(0,l) log(l-z) + A(0,0)| 

+ |A(i,l) log(l-z)+^(i,0)| Vl^+... (23) 

in terms of certain known coefficients which we denote here as A{m, n). For instance, A(0, 0) 
is the constant called K^'^'> in reference |15| . In terms of the conformal variable u) ([8]), this 
expansion reads (llOp 



n(2)(^) ^ 1,0) +|^(+)(o^2) log2(l-w) + f)(+)(0,l) log(l-a;) + f]W (0,0)1 

1 — uj y J 

+ |l7W(l,0) + 17W(l,l) log(l-L^)|(l-a;) + ... (2 



Obviously the two expansions (p3|) and (|24p are related, with 

-2 A(-i,0) 
4 yl(0,2) 

2 yl(0, 1) -8 ^(0,2) log 2 
)(0,0) = ^(0,0) + A(-i,0)-2^(0,l) log 2 + 4 yl(0, 2) log^ 2 



o(+)(- 


-1,0) 




(0,2) 




(0,1) 




(0,0) 




(1,0) 




(1,1) 



(1,0) =^(0,1)-- A(i,0)-4^(0,2) log2 + ^(i,l)log2 



(25) 



As we have discussed (jl3p . it is convenient to re-express the expansion (|24p in the t variable, 
(recall that a; = e~*), as 



n(2)(, 



j^(+)(-i,o) 



oj= e^) ~ 2) log2f+l^W(0,l) logt+JlW(0,0)+-l^(+V-l,0) 

i-i-O t 2 



+ |^J]W(-1, 0) - ^1^W(0, 1) + J^(+)(l, 0) + (l7(+)(l, 1) - 17(+)(0, 2)) logtj t + . . . . 



(26) 
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This means, through the Converse Mapping Theorem ()17ll8p . that the Mehin transform is 
given by 

M [n^^^e"*)] (s) = r(s) G+[s) 

_l^(+)(-l,0) 2r?(+)(0,2) J7(+)(0,l) f)(+)(0,0) + i 1,0) 
s — 1 s'^ s"^ s 

(0,2) (1,1) ^l^(+)(-l,0)-l l^(+)(0,l) + ^](+)(l,0) 

+ {s + iY + . + 1 ^^^^ 

where the last expression corresponds to its singular expansion. As one can see, the residues 
in this singular expansion are given by the coefficients of the threshold expansion (IMj) in a 
very specific way. 

Turning now to the OPE ([3]), we have [8] 

n(2)(z) ~ (5(0,2) log2(-4z) + 5(0,1) log(-4z) + 5(0,0)l 

+ |i?(-l,2) log2(-4z) + 5(-l,l) log(-4z) + i?(-l,0)l - 



z 



+ {B{-2,2,) log3(-4z) + 5(-2,2) log2(-4z) + 5(-2, 1) log(-4z) + 5(-2, 0) [> ^ + 



(28) 



(2'] (2) (2) 

where the constants 5(0, 0), 5( — 1, 0) and 5(— 2, 0) are called Hq , H)^ and ^2 in reference 
|15| . In fact, many more terms of the OPE are known [llj, but we stop at 0{z~'^) because it 
is enough for our illustrative purposes. 

In terms of the conformal variable oj ^ , one has 

n(2)(tj) ~ ^ |f)(-)(0,2)log2(l + w) + f)(-)(0,l) log(l + a;) + O(-)(0,0)| 
+ |f)(-)(l,0) +f)(-)(l,l) log(l + w)|(l+a;) 
+ |j7(-)(2,0) +J](")(2,1) log(l + tj) + r2(-)(2,2) log2(l + w)| (l + o;)^ + • • • , (29) 

where, again, the two expansions are related through 

O(-)(0,2) = 4 5(0,2) 

Jl(-)(0, 1) = -2 5(0, 1) - 16 5(0,2) log2 

(0, 0) = 5(0, 0) + 4 5(0, 1) log 2 + 16 5(0, 2) log^ 2 (30) 

and similar expressions for 0(-)(l, 0), (1, 1), (2, 0), (2, 1) and 0(-)(2,2) which, 
because of their length, we do not write out explicitly. 

Finally, in terms of the t variable (recall that for the OPE one has that oj = — e~*), this 
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expansion reads 



n(2) {uj = - e"*) I n^-^ (0, 2) iog2 1 + (0, 1) log t + o(-) (0, 0) 

+ {^J^(-)(o, 1) + ln(-\o,2) - ^i^(-)(i,o) - + f](-)(2,o) 

+ ^J7(-)(0,2)-^l^(-)(l,l) + f](-)(2,l) logt + f^(-)(2,2) log^tj + . (31) 

We emphasize that only even powers of t may appear in this expansion. Consequently, the 
Mellin transform (I17|18p becomes 



M 



n(2)(_e-*)l (.) = r(.) G^{s) 



l^(-)(0,2) 17(-)(0,1) , l}(-)(0,0) 

S'^ s 

fi(-)(2,2) 1 6 - 12 0(-)(2,l) -O(-)(0,2) 

(s + 2)3 +12 (s + 2)2 

1 O(-)(0,l) + 6 r?(-)(0,2) - 12 O(-)(l,0) - 12 0(-)(l,l) + 24 n(~\2,0) 
+ 24 JT2 ^^^^ 

where the last expression corresponds to its singular expansion. Notice how the poles are 
located only at negative even numbers (besides zero) in correspondence with the even powers 
in the t expansion. Again, the residues of the singular expansion are given by the coefficients 
of the OPE in (j29p in a very specific way. 

Following the discussion presented in section 2, the two singular expansions (I27p and (I32p 
can be obtained by realizing that the ^}{n) coefficients admit an expansion for large values of 
n of the form 



n^^{n) = ao,o + jai,o + ai,i lognl- + 02,0^ + O ( ^log^^ 

+ (-1)" |/3i,o + log^}^ + |/33,o + /33,i log"};^ 

+ |/55,o + /35,i logn + 135,2 log' r?j^ + o(Jj log^^ r?j , (33) 

where the powers of the logarithms log^^'^^ n correspond to the terms in the OPE and thresh- 
old expansion which have not been yet calculated. With this expression for O'^'^(n), it is 
immediate to compute the G±{s) functions introduced in (jl5ll6p . First we split the sum as 



n] n 



^(±)" n^^{n) n-' + ^(±)" [^{n) - O^^(n)] n"^ . (34) 



n=l 



n=l 



n=l 



Then, since some of the il(n) coefficients are known exactly, say for n = 1, ...,N*, one can 
actually approximate the second sum as 



N* 



G±(s) = ^(±)" n^'^{n)n-' + ^{±T ['^{n) - n^^{n)] n"^ + 5G±{N*,s) . (35) 



n=l 



n=l 
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where 

1 

5G±{N*,s) = ^ dt t^-i e{N\ ±e-*) , (36) 

and 

oo 

8{N\u:)= [n{n) - (n)] Lo'' . (37) 

n=N* + l 

In other words, up to the T{s) factor, 5G±{N* , s) in (j36p is the Mellin transform of the 
"error" functiorH i5(A^*,a;). The above expression for £(N*,uj) makes it clear that the error 
in the approximation will naturally become smaller as A^* gets larger, or as more terms in 
the expansion (pl) for are included. This is true over the whole physical region in the 
complex plane |a;| < 1. 

Using the expansion Q^^{n) in Eq. (|33p one can readily compute G-^- as 

G+(s) = ao,o C(s) + ai,o C{s + 1) - ai,i C'(s + 1) + "2,0 C(s + 2) 

+ /3i,o r/(s + 1) - r/'(s + 1) + /33,o r]{s + 3) - /33,i r/'(s + 3) 
+ /35,o ??(s + 5) - /35,i f?'(s + 5) + /35,2 r/"(s + 5) 

N* 

+ ^ [Jl(n)-f]^^(n)]n-" + (5G+(iV*,s) , (38) 

n=l 

and G_ as 

G_(s) = /3i,o C(s + 1) - C'(5 + 1) + /33,o C(5 + 3) - /33,i C'(5 + 3) 
+ /35,i C(s + 5) - /35,i C{s + 5) + /35,2 C"(s + 5) 
+ "0,0 ??(s) + ai,o ??(s + 1) - ai,i ??'(« + 1) + "2,0 ^/(•s + 2) 

N* 

+ ^(-1)" [0(n) -f)^^(n)] n~^ + (^G_(iV*,s) , (39) 

n=l 

in terms of the Riemann function and the Dirichlet rj function (recall Eqs. (j2ip ) and their 
corresponding first and second derivatives. Matching the two expressions (I27p and (I32p to the 
result in Eqs. ([38|) and (j39p one readily obtains the following results for the coefficients aij 
and Pi J in ([33]) : 

ao,o = 2 ^(-i,0) (- 3.44514) 

ai,o = -2 ^(0, 1) + 8-fEA{0, 2) + 8 log 2 A(0, 2) ( ~ -0.492936) 
ai,i = 8A(0,2) (=2.25) 
02,0 = ^(^1) (- 3.05433) 

"^An alternative expression is, of course, SG±{N* , s) = [f^Ci^) — n^'^(n)] n"", as one can see 

in (|35p . The important point is that this expression has to be understood by analytic continuation in s. 
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A,o = 2 B{0, 1) + 8 7i7;5(0, 2) + 16 log 2 5(0, 2) ( ~ 0.33723) 
/3i,i = 8 5(0,2) (- 0.211083) 

/33,o = -^(1, 1) + 6 5(-l, 2) - 47s5(-l, 2) - ^ 5(0, 2) - 8 log 2 5(-l, 2) ( ~ 0.183422) 
/33 1 = -4 5(-l, 2) ( ~ -0.620598) 

/35,o = 35(-2, 1) - 255(-2, 2) + 127^; 5(-2, 2) + 105 5(-2, 3) - 1507^^ 5(-2, 3) 

31 

+ 367I 5(-2, 3) - 6^2 B{-2, 3) - 5(-l, 1) + — 5(-l, 2) - A-fE 5(-l, 2) 

+ — 5(0, 2) + 24 log 2 5(-2, 2) - 300 log 2 5(-2, 3) + 1447^; log 2 5(-2, 3) 
15 

- 8 log 2 5(-l,2) + 1441og2 2 5(-2,3) (~ -4.30753) 
/35,i = 12 5(-2,2) - 150 5(-2,3) + 72 -fE 5(-2,3) + 1441og2 5(-2,3) - 4 5(-l,3) 
(~ -1.89016) 

/35 2 = 36 5(-2,3) (~ 1.38684) , (40) 

where (in parentheses) we also give their approximate numerical value (for 3 light flavors). 
We emphasize that these a and /3 coefficients in Eq. ([33|) are determined by the logarithms 
of the threshold expansion ([2]) and the OPE ([3]) in a totally unambiguous wajH 



0.(n) I 




5 10 15 20 25 30 



Figure 2: Comparison of the n^^{n) in Eqs. (|33p and (I40p with the first 30 exact coefficients 
n{n). 

Although the asymptotic expression (f33l) is, strictly speaking, only valid at large values 
of n, it is worth realizing that (j33p sets in very quickly, and it becomes a very efficient 
approximation even for relatively low values of n. For instance, for n > 2 — 3 the error is 
< 1%, and gets better for larger values of n. This could be compared to the value predicted 

^The coefScient ao.o is determined by the squared-root singularity in the threshold expansion (|23p . and is 
the only exception to this rule. This is all in agreement with the Converse Mapping Theorem [16| . 
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for n = 1, which has a much bigger error of the order of ~ 300%. To ihustrate this fact, we 
plot in Fig. [2] the exact value of the first 30 terms for Q{n) (obtained from (llip with the exact 
results for the first 30 C{n) from ref. O [ini El [7]) together with the result for the asymptotic 
expression (I33| l40l) . This high efficiency of the asymptotic expression (p3l) will be turned into 
a method to predict the original Taylor coefficients C(n) in ([1]) in the following. 

One can even get an idea of the O {n~^ log^^ n) error in the asymptotic expression (I33p by 
performing a fit to the tail of the values for C(n) for large n. Although the fit is very insen- 
sitive to the sign-alternating part of (I33p and, essentially, yields no information on the term 
0(n^''log^^ n), the equal-sign term 0{n~^ log^^ n) is numerically rather well approximated by 
a term of the form ~ +0.5 log^'^ n. Taking a 100% error on the value of the coefficient 
(but not on the sign) we estimate for the error: 

- i-rr o . (41) 

Although the second term in the previous expression is essentially unknown, the large power 
vJ in the denominator makes it essentially irrelevant for any reasonable value of ^2- This 
result will be useful for estimating the error of our calculations in what follows. 

3.1 Prediction for the constant K'^'^'> 

Given the result for G+(s) in (j38p it is straightforward to extract the residue at s = of 
the combination r(s)G+(s). According to Eqs. (llSMlSp . this residue determines the constant 
j4(0,0) in Eq. (j23p (which is called K^'^^ in ref. |15)). This prediction is interesting because it 
has not been possible to obtain the value for this constant by ordinary diagrammatic methods 
up to now. The result yields 

- ^ + + Y + ^i) "1.1 + Y«2,o - C'(2) a2,i - /3i,o log 2 + + 7^ log 2^ 

3C(3) / C(3)log2 3C(3)\ 15 / C(5)log2 15C-(5) \ 

P-ifl + + —r- P3,l - TT^ClSj /^5,0 + T7, + i« P5,l 



+ (CM _ _ 1^) ft, + £ [„(„) _ + eiN-^ I, , ,42) 

^ ^ n=l 

where 71 ~ —0.072816 is one of the Stieltjes's constants, appearing in the Laurent expansion 
of C(s) around s = 1. In general, one has that 



In = lim 



C(")(a;) - (-l)"■ 



(x-l) 



n+l 



(43) 



Although, in principle, the expression (I42p depends on the unknown term O (n~^ log^^ n) in 
Eq. (j4ip through £(N* , 1) in Eq. (pT]) . the numerical result for K^'^^ is largely insensitive to 
this, for any reasonable value of the leading power £2, due to the strong suppression of the 
n"^ factor in the denominator. Using Eq. (j4ip . our results in Eqs. (|40p and N* = 30 exact 
coefficients C(n) from [HIT], neglecting the second term O (n~^ log^^ n), one obtains 

= 3.783 i°:J]°^ . (44) 
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Ref. 




m 
m 

This work 


3.81 ±0.02 
3.71 ±0.03 

q 700 +0.004 
0.100 _o.OOO 



Table 1: Results for the threshold constant K^'^\ See text. 



This result is compared to two previous determinations, using Fade techniques, in Table [T]. 

We can now double check that the contribution from the O {n~'^ log^^ n) in Eq. (|4ip is 
indeed small, for any reasonable value of the power £2- In order to do this, let us note that 
analogously to the previous case, the residue at s = of the combination V[s)G-{s) also 
determines the value of the constant 5(0,0) (which was called Hq in ref. |15)). The result 
reads 



(2) 



"0,0 



lo ^ 2 ' 

log 2 ain ± ( 7£ log 2 — ) aij 



12 



2 /„2 
"2,0 ± 



TT" /7r2log2 C'(2) 



12 



± 



"2,1 



± ± ^ ± 71) ± C(3) /33,o - C'(3) /33,i ± C(5) /35,o - C'(5) /35,i ± C"(5) A 

N* 

+ ^(-1)" [n{n) - f]^^(n)] ± 8{N\ -1) . 



(45) 



n=l 



However, unlike the case of K^'^\ the result for h'^^ has been calculated[10j and, for 3 light 
flavors, it is given by 

= -0.58570 . (46) 

If we neglect the term O (n"'' log^^ n) in Eq. (j4ip . we can compute the value for the expression 
with the help of our results in Eqs. (|40p and the A^* = 30 exact coefficients C{n) {VL{n)). 



The result obtained is Hq [Eq. ^^ 



-0.5856 



+0.0001 
-0.0000 ' 



in agreement with (|46p . 



3.2 Prediction for C{n) coefRcients 



Expression (|12p gives us the exact dictionary to translate the values of the coefficients Q.{n) 
into those for C{n). Exact knowledge of a few C(n) (let us say for n = 1, A^*) together 
with the asymptotic expression J7^'^(n), allows us then to make a prediction for the rest of 
coefficients for n > A^*, as we will now describe. 

Splitting the sum in (|12p between 1 < p < N* and p > A^* and using the exact value of 
r2(p) for the ffi'st part and the asymptotic expression and r2^'^(p) (I33p for the second part, 
one obtains the following prediction C^^^'^{n): 



N* 



n{p) p 



^ r (1 ± n - p) r (1 + n + p) 

p=i 



Q'^^{p) p 



^ r (1 ± n -p) r (1 + n + p) 



(47) 
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Figure 3: Error of the approximation C^^^'^{n) in Eq. (j47p compared with the exact C(n), 
for different values of A^* = 2, N* = 3 and N* = 4, as explained in the text. Also shown is 
the direct asymptotic expression for C(n) obtained in Ref. [7] (labeled as "Asymp"). 



The values obtained for C^"^ {n) can now be compared, e.g., with the result of the direct 
asymptotic expression for C{n) obtained in ref. [7| (black curve in Fig. [3]). As one can see in 
this figure, already for A^* = 2 (i.e. only assuming that the two coefficients C(l) and C(2) are 
exactly known), one achieves a better approximation than that obtained directly in terms of 
the coefficients C(n) 0, with errors which are of the order of 6 per mil, at worst (blue curve). 
On this figure we also show the curves when A^* = 3 (and N* = 4), i.e. when 3 (respectively 
4) coefficients C{n) are exactly known. In these cases, the errors are even smaller (green 
(respectively red) curves). This shows that Eq. ()47|) together with the asymptotic expression 
(j331 [40|l is a good strategy to predict the value of the coefficients C{n). This fact will be 
relevant in the case of the function H^^). 

3.3 Reconstruction of n'^^) (2;) 
3.3.1 Analytic expression for Il^'^\uj) 

Starting from Eq. ([9]), we can again rearrange the series by singling out the first A^* coeffi- 
cients, again assumed to be exactly known, and write 

00 N* 

n(2)(a;) = '^n^s^n) + ^ [0(n) - O^^(n)] u:"" + 8{N\uj) , (48) 

n=l n=l 
^Except for n = 5, for which there seems to be a fortuitous cancelation. 
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where 

oo 

£{N*,oo)= Mn)-n^^{n)]u'' , (49) 

n=N*+l 

is the "error function" ah^eady defined in Eq. (|37p . The first term in Eq. (j48p can be exphcitly 
summed using the fact thai|f| 

oo , p 

Y'^u:- = {-iru(^Hs,u:), (50) 

Tl 

n=l 

where the polylogarithm may be defined in the w complex plane as 

Li(^''^) = T;7^ log^-^fi^ log(l-x^). (51) 



r(s — 1) Jq X \x 

Consequently, we find the following expression for n(2)(a;), 

n*^^) (a;) = aon ai o log(l - lo) - ai i Li^^^ (1, uj) + 02,0 Li(2, uj) 

I — UJ 

- /3i,o log(l +uj)- /3i,i Li(i)(l, -uj) + /33,o Li(3, -a;) - /33,i Li(i)(3, -u) 
+ 135,0 Li(5, -a;) - /^s.i Li^^) (5, -a;) + ^5,2 Li^'^ (5, -uj) 

N* 

+ J^[f)(n)-0^^(n)]a;" + £:(iV*,a;) , (52) 

n=l 

where the coefficients Oij and are given in (|40p . 

Neglecting the error function £(N*,u}) (an approximation which gets better as N* grows 
larger, and also as many more terms in the expansion (I33p are known), our expression for 
n(2)(a;) is given by a unique combination of polylogarithms (and derivatives thereof) plus 
a known polynomial of degree A^* in uj. This expression constitutes a resummation of the 
Taylor expansion ([9]) thanks to the information supplied by the threshold expansion and 
the OPE, as seen in Eqs. (I2]),([3|) and (|40|) . Notice that the knowledge of further terms in 
either expansion can only lead to the determination of further a's and /3's coefficients in 
^}^^{n), leaving unaltered the coefficients already determined. Therefore, the approximation 
is systematic and can be improved upon by either knowing more C{n) coefficients exactly 
(i.e. increasing N*) or by knowing more terms from the threshold expansion or the OPE 
(i.e. by knowing more a's and /3's in Q^^{n)). This should be clear from the fact that the 
expression for the error function £{N*,u;) in Eq. (j49p naturally becomes smaller in these two 
cases, for any \uj\ < 1 in the complex plane. In order to further illustrate this point, in the 
appendix we show the result of our resummation for the spectral function and compare it to 
the exact result, in the particular case of the function Il^^\z) in Eq. ([5]), where this result is 
fully known. 

3.3.2 Result for ImU^'^\z) 

Undoing the change of variables ([8]) , and using the expression for n(2)(a;) in one imme- 
diately obtains 

n..,.,.n<^.(i_^). ,53) 



^This function Li(s,a;) is nothing else than the usual polylogarithm Lis(aj). Here we write in this slightly 
unconventional way because we need to differentiate with respect to the variable s. We also use the notation 



Li'P'(s,tj) to signify ^Li(s,a; 
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Figure 4: Plot of the spectral function, Imll^^) i terms of the quark velocity v = 

y^l — 1/z. The dashed curves correspond to consecutive approximations from the threshold 
expansion [v — t- 0) and from the OPE — )• 1). 

as an representation for the function TI^^^ (z) in the whole z complex plane. In particular, the 
analytic continuation of the polylogarithms and the polynomial in co yield an approximation to 
its branch cut and the corresponding spectral function. In principle, we could easily give the 
explicitly analytic result for this spectral function from the expression (|521 [53]) but we refrain 
from doing so because the result is not particularly illuminating. Alternatively, we show in 
Fig. mthe result for this spectral function as a function of the quark velocity v = y^l — 1/z 
for N* = 3, together with the corresponding approximations at threshold {v — )• 0) and the 
OPE (v — )• 1) (dotted curves). In this plot, the error band is estimated with the help of the 
error function £{N* = 3, w) in (f52|) and the educated guess (|4ip . Even for A^* = 3 the errors 
are so small that they cannot be seen on the plot, becoming even smaller at larger values of 
N* . This plot may be compared to the one in Fig. 3 of ref. [15j. One could of course use all 
the known C{n) coefficients, i.e. take A^* = 30 in our expressions. However, at the current 
level of precision, this is not yet necessary |15) . 

4 The case of n(3)(g2) 

We can now apply the same method exposed in the previous section for the H^"^^ function, to 
treat the function H^^^q^). For later convenience, let us first reproduce here the coefficients of 
the Taylor, OPE and threshold expansions listed in ref. |13] for the vector vacuum polarization 
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function at ©(af) (for 3 light flavors), 
U'-^\z) = 6.95646 z + 7.24783 + 7.31859 + 0{z^] 



2.63641 25.2331 + 7.75157 log(l - z) 

z^l 1 — Z 



+ K(3) 

11.0654 log(l -z) + 1.42833 log2(l - z) - 0.421875 log^(l - 



+o(VT^ iog^(i-2) 



- 0.0698778 log(-4z) + 0.121085 log2(-4z) - 0.0366469 log^{-iz) 



2— > — OO 
1 



+ - -3.75665 log(-4z) + 2.11728 log2(-4z) - 0.318916 log^(-4z)| 



^ ^H^ - 5.13014 log(-4z) + 0.31818 log2(-4z) + 0.401495 log^(-4z) - 0.0786514 log''(-4 

+0 (z-^ log^'(-42)' 



Using the coefficients of these expansions, a straightforward rerun of the method leads to 
the following result for the asymptotic expression of the corresponding coefficients ^}^^{n) in 
this case: 



Q^^{n) = a_i,o n + \ ao,o + ao,i log n [> + <j ai,o + ai,i log n + ai,2 log^ n\— + O [ log*^ n 



1 



1 



+ (-ir 



/3i,o + log n + /3i,2 log2 n I i + | /33,o + (^3,1 log n + ^3^2 log^ n I ^ 

n I \ 



+ \ /35,o + /35,i log n + /35,2 log2 n + /35,3 log=^ n\^ + 0[^ log^^ ^ 



(55) 



witl 



a_i,o =i 10.5456 
ao,i ~ 31.0063 
ao,o ^ -11.0769 
ai^o — 36.3318 
ai,i ~ 37.1514 
ai2 ~ 10.125 



f /3i,o =i -0.181866 
~ -2.4852 
/3i,2 =i -0.879515 
/33,o - -10.4385 
/33,i ~ -4.7750 

L /33,2 =^ 3.82702 



r /35,o - -70.9277 
/35,i ~ 56.3093 
/35,2 =i 20.9951 

I /35,3 - -7.55063 



(56) 



This expression for il^'^(n) can now be compared with the exact numerical result for 
0(n) which can be obtained from the exact value of the ffist 3 coefficients C(n), n = 1, 2, 3, 
calculated in ref . |2H [22] [23l [23] , and our previous formula (jlip . The result is shown in Fig. 
[5j As one can see, the asymptotic form (|55| [56|) sets in very quickly, as it also happened in 
the previous case of the function II^^^ (see Fig. [2]), reproducing the exact values for $1(2) and 
ri(3) within only a few per cent error. Notice that the error for ^2(1) would be much larger 
~ 300%, as it is clearly seen in Fig. [5j Based on this success for the "prediction" of $7(2) and 



"Here we only quote the numerical results, for simplicity. 
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n^^(n) , n(n) 



150 



100 



50 



o 



10 



Figure 5: Comparison between the asymptotic expression n^^{n) in Eqs. ([55]) and ([56]) (blue 
crosses) and the exact value for the first three coefficients ^}{n), n= 1,2,3 (red circles). 



0(3), we can make a very rough guess for the error of the asymptotic expression (|55p in the 
form: 



[VL{n) - VL^^n)] 



n>Af*(-2-3) 



zb 15 



log n 



± (-1)" O 



log^^ n 



n' 



(57) 



At present we know nothing about the sign-alternating term ~ log^^ n in the expression 
above although, fortunately, one is protected by the high power in the denominator, i.e. n^, 
for any reasonable value of the ^2 power. This also happened in the case of the function H^^^ 
(see Eq. (I4ip ). Therefore, the error coming from the first term in (|57ll and, moreover, its 
slow fall-off with n, i.e. ~ n^^, will now be the main source of error in all our estimates. 
It is important to realize that the error (I57p can only be reliably reduced once the terms 
of C(n~^ log^^ n) in Eq. (j55p are computed, which in turn requires the calculation of the 
terms 0{^/l - z log^(l - z)) , (£ / 0) in the threshold expansion (the second of the Eqs. 
(j54p ). Regretfully, this information is not available at present. Alternatively, one could also 
compute more C{n) coefficients from the Taylor expansion, in order to increase the value of 
A'"*. However, this second possibility does not seem so promising due to the rather slow ~ n"^ 
fall-off with n shown by the non sign-alternating part of the expression in Eqs. (j55p and (|57p . 
At any rate, using the expression for il^'^(n) in ()55|56p . the result for the constant K^^^ 
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in Eq. (j54|n . akin to the constant K^'^^ in the section 3.1, now becomes 



"-1,0 



ao,o log 27r , /^^^ , 7s , ^ . f 1% ^^7i? _ , , , 



2C(3) \ 



4 2 ^2 

^ + ^ + ^ + ^3 + '^^Em ) ai,3 - /3i,o log 2 + 



log2 2 



+ log 2 /3i,i 



+ I - + 27i? log^ 2 + 271 log 2^ /3i,2 + f 7E log^ 2 - + 871 log^ 2 + 872 log 2 ) /3i,3 



8C(3) ^ ^ / C(3)log2 ^ 303) ^| ^ ^ |^ C(3)log^2 ^ C^(3)log2 ^ 3C"(3) 



+ 



4 ' V 4 4 y ' V 4 

V(3)log^2 _ 3C^(3)log^2 ^ 3^(3) log 2 ^ SC^^^^ 

,.2 



'3,2 



/33,3 



15C(5) 



\ 4 4 4 4 / ' ' 16 

/3C(5)log2 15a5)\ , ^ C(5)log^2 , C^(5)log2 1505) \ 

V 48 + V 16 ^ 8 i^j 

/C(5)log^2 3C'(5)log2 2 3C"(5)log2 15C(^H5)\ 



/35,0 
/35,2 



16 16 16 16 / 

'C(5)log^2 C'(5)log^2 3C"(5)log2 2 C,^^\b)\og2 14C(^)(5)\ 



16 



+ 



16 



5,4 



(58) 



n.=l 



where 7^ are the Stieltjes's constants defined in Eq. (j43j) . Using again the 3 exactly known 
coefficients Q.{n) (n = 1,2,3) from |13| (i.e. N* = 3) and our estimate for the error in the 
coefficients (|57p to compute the error function £{N* , 1) defined in ([37]), we obtain: 



K^^^ = 18 ± 87 



(59) 



Although our central value almost coincides with that in ref. |13j . our error is ~ 8 times 
larger than that quoted in this reference and, in particular, cannot exclude a negative value 
for as it was found in Ref. |15) . The situation is summarized in Table [2l It is now 

immediate to assess the impact of a reduction of the systematic error ()57p . For instance, if 
this error were reduced to ~ ±15 n"^ log n one would get an error in 

if (3)^ Sk(^) ^ ±10, 

whereas if the error could be pushed down to ~ ±15 n"^, the error would be ^ ±0.6; 

a dramatic improvement. 



Ref. 




m 
m 

This work 


17± 11 
-10.09 ± 11 

18 ±87 



Table 2: Results for the threshold constant K^^\ See text. 



■'This constant is called K^o ''" in Eq. (19) of ref. [13 
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n 


(167rV3)CP'''=^(n) for iV* = 3 


Error 


C(n) from [151 


Error 


1 


366.174 





366.175 





2 


381.510 





381.509 





3 


385.235 





385.233 





4 


382.7 


0.5 


383.073 


0.011 


5 


378.0 


1.2 


378.688 


0.032 


6 


372.5 


1.8 


373.536 


0.061 


7 


367.0 


2.3 


368.23 


0.09 


8 


361.5 


2.7 


363.03 


0.13 


9 


356.4 


3.1 


358.06 


0.17 


10 


351.6 


3.4 


353.35 


0.2 



Table 3: Table of comparison between the C^^ (n) from this work and the C(n) in [13j. 



(3) 

As a cross-check of our method, we can also calculate the value of the first constant Hq ' 



appearing in the OPE [15] . We obtain 



iff 



—r IT- H ^ «o,i - log 2 ai n + h 7s log 2 ai i 



2 

+ ( 271 log 2 - + -iE log^ 2^ ai,2 + (^371 log^ 2 + 372 log 2 - + 7^ log^ 2 ) ai,; 



3 

+ (1^ + ^ + ^+73 + 27sC(3)) /3i,3 + C(3) /33,o - C'(3) /33,i + C"(3) /33,2 

- C^'^(3) /33,3 + C(5) /35,o - C'(5) /35,1 + C"(5) /35,2 - C^'^(5) /35,3 + C^'^(5) /35,4 
N' 

+ Y,[n{n)-n^^{n)\{-lY + £{N\-l) . (60) 

n=l 

Using the values obtained in (|56]l . our estimate (l57|) and A^* = 3, we obtain W^' = -7.1 ±1.2 
while the exact value is h'^^ = —6.1717 [13j. The reason why this estimate is much better 
than (|59p is because K^^^ is more sensitive to the non sign-alternating component of the error 
(j57p . which is large. Further constants, such as H^'' and -^2'^^ in (|54p . are more sensitive to 
this component of the error, due to the extra powers of n in the sums, resulting in a larger 
total error and we refrain from giving their result. 

Concerning our prediction for the coefficients C^^'^'^{n) , we can use the exact value for 
the ffist three C(n), n = 1, 2, 3, again from [HI EU [23l [2l] , and our result in (gTI with the 
expression for Q^^{n) from (|55l) and the values (|56l) to obtain the results listed on Table 
[31 i We find agreement with the numbers quoted in [13] . although our central values are 
systematically lower and the errors are larger. 

Finally, we can also give the expression for II^^) obtained by summing the coefficients 
(|55|56p . in complete analogy to what we did in Eq. (I48p . The result is 



^There is factor IGtt'^/S of difference between our normalization and that in Ref. [13) . 
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0.2 0.4 0.6 0.8 1.0 

Figure 6: Result for the spectral function of (|6ip in terms of the quark velocity, v, together 
with the threshold approximation {v — t- 0) (in green) and the OPE (v — )• 1) (in red). 



n(3) (cj) = a _i,o " ^9 + ao,o ao,i Li^^^ (0, uj) - ai,o log(l - u) 

(^1 — UJJ'^ 1 — UJ 

-ai,iLiW(l,w) + ai,2Li(2)(l,a;) -/3i,olog(l + tj) -/3i,iLiW(l,-L^) 
+ /3i,2 Li(2) (1, -a;) + /33^o Li(3, -w) + /33,2 Li(2) (3, -w) + /S^^o Li(5, -a;) 
- /35,i Li(i) (5, -w) + /35,2 Li(2) (5, -a;) - /35,3 Li^^) (5, -^^) 

AT* 

+ J2[^(.n)-n^^{n)]uj'' + £{N*,uj) . (61) 

n=l 

Again, we plot in Fig. [6]the result for the spectral function which comes from the imaginary 
part of (j6ip . This figure may be compared to the akin figures in ref. |15| (Fig. 5) and ref. 
[T3] (Fig. 2, first panel). 

5 Conclusions and Outlook 

Given an analytic function with a cut located at 1 < z < cx), we have presented a method for 
reconstructing the full function from its Taylor expansion around the origin, i.e. z — )• 0, Eq. 
([T|), its threshold expansion for z — )• 1, Eq. ([2]), and its OPE series for z — )• oo, Eq. ([3]). We 
have exemplified this method with the case of the perturbative vacuum polarization function 
of a quark, but we think that it may be useful in other contexts as well. 
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The method uses the conformal variable u, Eq. ([8]), in terms of which the whole cut 
complex plane is mapped onto a disk of unit radius (with the cut running at the boundary, 
see Fig. [1]), and consists in realizing that, in this variable, the Taylor coefficients Q{n), Eq. 
([9]), admit an asymptotic expansion in inverse powers of n, for large n, with a sign-alternating 
component and an equal-sign component. The coefficients of the equal-sign component, which 
we call ap^k in Eq. (I22p . are essentially determined by the logarithms of the threshold ex- 
pansion, whereas the coefficients of the sign-alternating component, which we call /3p ^ in 
Eq. (I22p . are determined by the logarithms of the OPE. Since logarithms are controlled by 
the Renormalization Group, it is not inconceivable that one may also link this asymptotic 
behavior in n to the Renormalization Group in the future. 

We have found that the precise dictionary, Eqs. (I17|18p . is most easily established by 
means of the use of the Mellin transform and the Converse Mapping Theorem exposed in 
Ref. [16] and references therein. Once an asymptotic expression for the coefficients Cl^^(n) is 
obtained, the power series can be explicitly summed using this asymptotic expression starting 
from a certain A'^* onwards. Below this A^*, the exact values for the coefficients must be 
used. The precision of the result of course depends on how precise the value of ^^^{n) is for 
n = N* . Interestingly, we have found that the asymptotic series Q^^{n) sets in very quickly, 
giving quite accurate results already for n > 2, 3. 

The resummed function in terms of the conformal variable oj can then be expressed as a 
combination of polylogarithms and their derivatives plus a polynomial, with known coefficients 
(see, e.g. Eqs. (I52p . (l6T]) : and the Appendix). The method is analytic and it is possible to 
know how our results depend on our ignorance due to uncalculated terms in the Taylor, OPE 
and threshold expansion; which has given us a way to estimate the associated systematic error 

m- 

We have made some numerical comparisons with the results obtained via the Pade method 
[Ml [131 [H] , ™d our results agree with them, although our errors are usually larger. Never- 
theless, we think it would be very interesting to be able to make a more thorough comparison 
including, in particular, the result for the full spectral function. Furthermore, we plan to use 
our method to study other channels as well [26| . 

Regarding the 0{a^) vacuum polarization function, II^^^ we find that the calculation of 
all the terms 0{\/l — z log^(l — z)) from the threshold expansion for ^ 7^ 0, or at least the one 
with the highest power i, would result in a substantial reduction of the associated theoretical 
error in all the predictions for this function. 
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APPENDIX 

In this Appendix we will use a known function to further illustrate our method. For this 
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5 10 15 

Figure 7: Relative error between the asymptotic form il^'^(n) 
values for Q(n) as a function of n. 
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in Eq. (l62ll and the exact 



purpose, we will choose the 0{a^) result, 11(^^(2:), which reads explicitly 



n(o)(z) 



167r2 



20 4 



4(l-z)(l + 2 z) 



3 z 



G{z) 



where 



and 



Giz) 



2 u logu 



yr^T/i-i 



yr^T7i+i 

Using the method we describe in the text, the OPE and threshold expansions for this 
function determine that the asymptotic expansion of the associated coefficients 0(n) in the 
Taylor expansion of n(o) (w) in the conformal variable uj are given by 



n 



AS, 



n] 



1 



9 



1 

27r2 n 327r2 



— + 



(62) 



Amusingly, there are only sign-alternating coefficients, i.e. but no atj coefhcient. This 
feature is in agreement with the fact that the threshold expansion of Yl^'^^z) does not have 
any log(l — z) term, so all the terms in ^l''^^{n) are given by the \og{—z) terms of the OPE. 

As also happens for Jl^'^\z) and Ii^^\z), the asymptotic expansion Vl^^{n) sets in very 
quickly. In Fig. [7] we compare the exact values for the coefficients 0(n) and this asymptotic 
expression. As one can see, already for n = 2 the error is only a few percent, becoming smaller 
for larger n. 

As to the approximation for the full function n(^)(z), we obtain, in terms of the oj variable, 
the following expression: 
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imn«"(^ 



0.001 




Figure 8: Relative error in the spectral function for the approximation H^'^) and the 

exact value H^^^ (z) as a function of the quark velocity v = y^l — 1/z. The three curves shown 
correspond to A^* = 2 (blue), A^* = 4 (red) and A^* = 6 (green) 



N* 



n(0)(u;)Upp. = J]f]^^(n)^" + ^[J7(n)-J7^^( 

n=l n=l 

= ^ log(l + ^) + ^ Li(5, -w) + Yl [^(^) - 



N* 



n)\ oj 



n=\ 



which, as usual, depends on how many coefficients n{n) we assume to be exactly known 
(N*). In Fig. [8] we show how this expression approximates the exact function over the whole 
complex plane by comparing the spectral function obtained from this approximation to the 
exact spectral function of the function Il^^\z), for different values of Ar*oAs we can see on 
this plot, even for A^* = 2 (i.e. assuming only that the first two Q{n) coefficients are exactly 
known) we reproduce the exact result with a relative error at the per mil level (or smaller), 
over the whole range in velocity < v < 1. For larger values of A^* the errors are, of course, 
even smaller. 



^Due to the dispersion relation ([6}, it is clear that approximating the spectral function is equivalent to 
approximating the full function in the whole complex plane. 
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